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We develop here a general formalism for multi-orbital Mott systems which can be used to under- 
stand dynamical and static spectral weight transfer. We find that the spectral weight transferred 
from the high energy scales is greatly increased as a result of the multi-orbital structure. As a 
consequence certainly dynamically generated symmetries obtain at lower values of doping than in 
the single-band Hubbard model. For example, in the atomic limit for filling less than one elec- 
tron per site, the particle-hole symmetric condition in the lower band shifts from the one-band 
result of x = 1/3 to x = l/(2n + 1), where n is the number of orbitals with an unpaired spin. 
Transport properties computed from effective low-energy theories which forbid double occupancy 
of bare electrons, such as the multi-orbital t-J generalization, should all be be sensitive to this 
particle-hole symmetric condition. Away from the atomic limit, the dynamical contributions in- 
crease the transferred spectral weight. Consequently, any phenomena which are sensitive to an 
emergent particle-hole symmetry should obtain at x < (l/(2n + 1). 

PACS numbers: 71.27. +a, 74. 70.-b 
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I. INTRODUCTION 

The parent of the copper-oxide superconductors is a 
single-orbital Mott insulator in which one unpaired d- 
electron is responsible for the essential physics. By con- 
trast, the iron pnictides are multi-orbital bad metallic 
systems in which several d-electrons play an active role. 
Given that both systems superconduct, it is highly sug- 
gestive that combining the key ingredients of the two 
might result in a higher T c material. One of the ubiqui- 
tous features of the cuprates is dynamical spectral weight 
transfer— ~— arising fundamentally from the lack of rigid- 
ity of the Hubbard bands in a Mott system. Such spec- 
tral weight rearrangments are well documented in op- 
tical conductivity measurements-^ and oxygen K-edge 
x-ray absorption 3 . Spectral weight rearangements in 
the cuprates are well described^ by the single-band 
Hubbard model. Classic Mott systems such as VO2 
and NiO contain several d-orbitals that are singly oc- 
cupied. At present, although there are some studies on 
Mott transitions in multiorbital systems using dynami- 
cal mean-field theory (DMFT)^"— , there is no formal- 
ism that describes the spectral weight rearrangements 
in such systems under doping nor any experiments that 
observe this effect. In particular, since DMFT neglects 
the momentum-dependence in the single particle self- 
energy, the essential physics of hopping-induced spectral 
weight transfer under the introduction of dopings can 
not be captured correctly. Such momentum dependence 
is particularly important in the treameni^ of mixed- 
valent systems in which inter-site coherences dominate 
the physics. It is precisely the theoretical description 
of spectral weight rearrangement in multi-orbital Mott 
systems that we tackle here. We develop the strong cou- 
pling perturbative analysis for multi-orbital systems and 
derive expressions for the spectral weight shifts in the 
lower band as a function of doping. The Hund's couling, 
a new energy scale absent in the single band Hubbard 



model, is included explicitly in this study. It has been 
pointed out that Hund's coupling play a crucial role in 
determining the correlation effects^. As a result, it is 
also important to undersatnd how Hund's coupling af- 
fects the spectral properties which is also investigated 
in this paper. The key result we obtain is that multi- 
orbital physics enhances the total spectral weight trans- 
ferred across the Mott gap. It is our hope that this study 
will motivate experiments on doping multi-orbital Mott 
insulators. 

Indeed, the low energy physics of strongly correlated 
systems remains one of the most challenging problems. 
The difficulty in strong coupling is uncloaking a set of 
excitations which get rid of the interactions. The pres- 
ence of spectral weight transfer demonstrates that such 
entities, should they exist, are not straightforwardly re- 
lated to those in the UV-complete theory. For example, 
it has been known since 1967— in the context of the 1- 
band Hubbard model that the total spectral intensity 
of the lower band exceeds the total number of electron 
states in the band. Hence, the additional degree of free- 
dom in the lower band must be orthogonal to an electron. 
The culprit that accounts for the excess spectral weight 
in the lower band is dynamical spectral weight transfer 
(DSWT)^ii2. DSWT emerges from mixing between dou- 
bly occupied sites in the upper band with singly occupied 
ones in the lower band. The fact that this effect can ac- 
count for the orthogonality for electron addition at low 
energies implies that it provides a natural mechanism 
to explain the vanishing of the quasiparticle weight, Z, 
seen in angle-resolved photoemissionM (ARPES) in the 
underdoped cuprates. Consequently, it is important to 
catalogue precisely how this effect is manifest in multi- 
orbital Mott systems. We demonstrate here that this ef- 
fect is even more pronounced because the multi-orbitals 
enhace the phase space for dynamical mixing between 
the electronic states below and above the Mott gap. 

This paper is organized as follows. In Section II, we 
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consider a general multi-orbital Hubbard model in the 
presence of Hund's coupling and analyze the spectral 
weight transfer for various electron fillings in the atomic 
limit. Due to the existence of Hund's coupling and the 
paired hopping with the energy scales J, we show that 
the distribution of the spectral weights is enriched. More- 
over, the non-zero J also complicates the calculation of 
DSWT when the hoppings are introduced. This is an- 
alyzed in Section III. The full details of the derivations 
are placed in the Appendix. 



II. SPECTRAL WEIGHT TRANSFER IN 
ATOMIC LIMIT 

The Hamiltonian describing a multi-orbital Mott sys- 
tem consists of two parts, 

H = H K +H X , (1) 

with Hk the non-interacting tight-binding Hamiltonian 

h k = - > ; t ab 4 aa c Jba , (2) 



E 

(tj)aba 



and c\ aa creates an electron with spin a in the orbtial 
a on site i. Due to the existence of the multi-orbital 
structure, electrons on the same site experience differ- 
ent types of interactions, including the intra- and inter- 
orbital Coulomb repulsive interactions, Hund's coupling, 
and the paired hopping between different orbitals. The 
resulting on-site interaction for the multi-orbital system 
reads 



Hi = J^Un 



,7 



i.b>a 



^ 2 J S ia ■ Sa + J (PiaPib + h.c.) , (3) 



i,b>a 



where 
lJ 



n ia\ + n ia\,, Si a — 



state from the high energy sector, thereby leading to the 
1 — x residue. The empty state reappears at low energies 
as part of the addition spectrum at low energies. The 
fact that the empty part of the spectrum grows as 2x 
as opposed to x in a band insulator reflects this static 
transfer of spectral weight. 

To begin our study of a multi-orbital system, we start 
with the two-orbital case since a variety of transition 
metal oxides with orthorhombic crystal structure have 
degenerate quasi- ID d xz and d yz bands which dominate 
in the low energy physics. In this case, the Mott physics 
will be prominent near integer fillings of n = 1, 2, 3. One 
interesting question is how is the spectral weight redis- 
tributed when such systems are doped with holes or elec- 
trons near these integer fillings. In the atomic limit, the 
spectral weights can be obtained exactly by evaluating 
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where 



E 



E 



JV±1 



E% =F ir, 



(5) 



(6) 



is the electron removal (addition) Green function on site 
i, and and {V'ml are the eigenenergies and eigen- 

vectors of Hi with two orbitals listed in Table [U Because 
a non-zero J lifts the degeneracy of two-electron states 
on a single site, the spectral weight also splits into sev- 
eral parts separated by the energy scales of J whenever 
two-electron states are involved. For example, if J = 
and n = 1 the spectral intensity should be 3 for the 'up- 
per Hubbard bands' and 1 for the 'lower Hubbard band' 
with a gap of U between them. As J is turned on, using 
Eq. [5] and Table [H we have that 

A+H = ls{u - /x - U' + J) + l -6{u - n - U' - J) 
+ ^5(lj - fi-U+J) + ^S(uj -fi-U-J), 



2 c ian°nv c ia.v, Via = CiafCiai- It should be noted that A i (w) = S(u - fi). 



(7) 



U' = U — 2J holds for systems with both orbital degen- 
eracy and cubic symmetry. This relation will be used in 
this paper unless stated otherwise. 



A. Two-orbital Mott system 

Before we present the result for the 2-orbital system, 
let us recount the atomic-limit result for the single-orbital 
Mott insulator. In this limit, the exact single-particle 
retarded Green function is 



Gf(w) 



1 



1 



u 

2 



U_ 

2 



(4) 



The two poles, split by U, carry spectral weight of 1 + 
x and 1 — x respectively. Each hole subtracts a single 



It can be seen that the upper Hubbard bands split into 
three parts with energies of U — 3 J, U — J, and U + J 
(V = U — 2 J is used) with different spectral intensities 
as shown in Fig. (Ua). As the system is doped with ei- 
ther holes or electrons, the redistribution of the spectral 
weights occurs at all energy scales in a manner similar 
to that in single-band Hubbard models. Moreover, the 
amount of spectral weight transferred from the high en- 
ergy down to low energy is even larger in multi-orbital 
Mott systems. The reason is that the multi-orbital struc- 
ture provides more ways to add or remove electrons. In 
the following, we will analyze case-by-case the spectral 
weight distributions near different interger fillings. 

The case of n — 1 — x is the most prominent example to 
illustrate the effect of multi-orbital structure on spectral 
weight transfer. For an Af— site system, there are Af(l — 
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TABLE I: A complete list of the eigenstates and eigenenergies 
of Hi given in Eq. [3] for two orbitals on site i. 


No. of electrons (N) 


Eigenstate 


Eigenenergy 


1 


Any 1-electron state 





2 


|lt,2t> , |H,2|) 
^(|lt,2|) + |H,2T» 
^(|lt,2|}-|H,2T» 

^(|it,H)-|2t,2;» 

^(|lt,H) + |2t,2;» 


U'-J 
U'-J 
U' + J 
U-J 
U + J 


3 


Any 3-electron state 


U + 2U' -J 


4 


|lt,H,2t,2|) 


2U + AU' -2J 



x) sites with one electron and Afx sites empty in the 
ground state. For sites occupied by one electron, Af(uj) 
are the same as in Eq. [JJ For unoccupied sites, we find 
that 



A+(cj) = 4tf(w-/i), 
At(uj) = 0. 



(8) 



Therefore the intensity of the lower energy bands near 
the chemical potential is simply 

m h (l-x) = duA7(u>)+ fduA+iu)} 

; J — OO J U 



1 - x + 4% = 1 + 3x, 



(9) 



where A is a cutoff energy scale demarcating the division 
between IR and UV scales. One can realize Eq. [9] by 
counting the total number of ways to remove and add 
electrons to the system without resulting in an energy 
change. While there is one way to remove an electron 
from an occupied site contributing 1 — x to m (1 — x), 
there are four ways to add an electron to an empty site in 
a two-orbital system with spin 1/2 leading to the remain- 
ing contribution of 4x. This is precisely the same argu- 
ment used to understand the spectral properties in the 
single-band Hubbard model leading to the well-known re- 
sult (see Eq. (j4|)) that the total weight of the lower band 
is rriLHB = 1 — x + 2x = 1 + a4^, where x is the number 
of holes. In fact, for n = 1 — x, the spectral weight can be 
obtained straightforwardly based on the same argument 
for n -orbital systems with spin 1/2 in the atomic limit, 



m' l (l 



1 



+ 2n a x = 1 + (2n - l)x. (10) 



Clearly, the multi-orbital structure amplifies the spectral 
weight transfer from high energies down to the chemical 
potential even in the atomic limit. 



B. n -orbital Mott system 

The generalization of the spectral weight transfer for 
n > 1 is complicated by Hund's coupling and paired hop- 
ping terms which disperse the spectral weights at energy 
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FIG. 1: The schematic illustration of spectral weight transfer 
in a two-orbital Mott system for (a) n = 1, (b) n — 1 — x, and 
(c) n = 1 + x. U' has been replaced by U — 2J in the plot. 
The y-axis represents the magnitude of the spectral weight 
and the x-axis, the energy. Although the bands are strictly 
S— function peaks in the atomic limit, they are represented as 
broadened here entirely for aesthetic purposes. The dashed 
line represents the location of the chemical potential. 



scales on the order of J in the lower and upper Hubbard 
bands. Nevertheless, since the Mott insulating phase at 
integer filling requires a well-defined Mott gap Anjott de- 
marcating the lower and upper energy bands, it is ex- 
pected that the spread of the spectral weight due to J 
near the chemical potential occurs within the Mott gap. 
In other words, we can always find proper cutoff energy 
scales fx ± E ± such that E ± ~ O(J) < A Mo tt- Then for 
hole-doping near integer fillings in an 7i -orbital system, 
the intensity of the bands near the chemical potential can 
be evaluated using 



m [rii — x) 



duA+(u)] 



= TLi{\ — x) + (2n — TLi + l)x 

= m + (2n - 2rii + \)x. (11) 

The above result can be understood as follows. The 
ground state with filling n = n.j — x has a total frac- 
tion of (1 — x ) sites occupied (per site) by rii electrons 
and a total fraction of x sites occupied by rii — 1 electrons. 
The processes with energies on the order of J away from 
the chemical potential correspond to removing electrons 
from sites occupied by nt electrons and adding electrons 
into sites occupied by rii ~ 1 electrons. Consequently, 
there are rii ways for electron removal and (2n — m + 1) 
ways for electron addition. Eq. [IT] just accounts for the 
sum of these ways to add and remove electrons discussed 
above. 
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FIG. 2: Schematic illustration of spectral weight transfer in 
a two-orbital Mott system for (a) n = 2, (b) n = 2 — x, and 
(c) n = 2 + x. 17' has been replaced by U — 2 J in the plot. 



III. DYNAMICAL PECTRAL WEIGHT 
TRANSFER 



In this section, we proceed to compute the dynamical 
spectral weight by treating Hk as a small perturbation 
in the strong-coupling limit (t a b « U, U'). The method 
used here is a generalization of the unitary transforma- 
tion method^ i 13 i 15 i 16 which has been successfully applied 
to study DSWT in the single-band Hubbard model. The 
formalism given below is general regardless of the mag- 
nitude of V = U — 2 J. Therefore, we will write U' as an 
independent parameter from now on. 



A. General formalism of the perturbation theory 

The central idea of this method is to perform a unitary 
transformation 

H = e s He- s , H = H K + Hi , H = H K + Hi (14) 

with the condition 



Likewise, for electron-doping, n = n% + x, the ground 
state is composed of Afx of total sites occupied by rn + 1 
electrons and Af(l — x) of total sites occupied by m elec- 
trons. As a result, there are + 1 ways to remove elec- 
trons from sites occupied by r^ + l electrons and (2n G — rij) 
ways to add electrons into sites occupied by n, electrons 
with energies of the order of J away from the chemical 
potential. This yields 



m e (rii + x) 



I 



li- E- 



cL,A+(u,)] 

(m + l)x + (2n - rii){l - x) 
2n Q — rii + (1 + 2rii — 2n )x. 



(12) 



The spectral weights evaluated from Eq. [5] and Table U 
for fillings of n — l±x and n = 2±x in a two-orbital Mott 
system are illustrated in Figs. Q] and [2] Proper choices 
of cutoff energies (E~,E + ) to account for the spectral 
intensity near the chemical potential m e,h (n = m±x) are 
(-co, 0+) for n = 1 - x, (0" , 4 J) for n = l+x, (0" , 4 J) 
for n — 2 — x, and (4 J, + ) for n = 2 — x respectively. It is 
interesting to mention that this generalized intensity of 
the bands near the chemical potential is symmetric with 
respect to 'half-filling' (rij = n ), i.e., 



m h (rii — x) — m e (2n — m + x). 



(13) 



[H, Hi] = 0. 



(15) 



In other words, the goal is to find 'dressed' fermion opera- 
tors Ci aa - such that the original Hamiltonian H shares the 
same eigen basis with the transformed interaction Hamil- 
tonian Hi. Since Hi conserves the number of multiply 
occupied sites in terms of these dressed fermionic degrees 
of freedom, so does the original Hamiltonian H . As a re- 
sult, the ground state of H can be expressed in terms of 
the number representation of the dressed fermionic de- 
grees of freedom. In general, for any operator O, we 
define O such that 



= 0(c), 6 = 0(5), 



(16) 



simply by replacing the Fermi operators Ci aa with the 
transformed fermions Cj aiT . 
Expanding Eq. [T4l yields 



H = H+[S, H] + ^[S, [S, H\] 



(17) 



S is designed to satisfy Eq. [15] order by order in pertur- 
bation theory 



S = Si + s 2 + 



(18) 



which are all written in terms of dressed fermionic degrees 
of freedom. The first-order Si solves 



[[Si, - [-H K ,Hi] 



(19) 



A naive solution could be [Si, Hi] = —Hk, but further 
simplifications can be made as explained below. In prin- 
ciple one can always classify Hk according to the change 
of the interaction energy for the corresponding hopping 
processes. Hopping processes conserving the interaction 
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energy, which is denoted as H K , commute with Hi and 
thus produce no important effects in this perturbation 
expansion. As a result, the suitable choice for Si should 
be 



[S 1 ,H I ] = -(H K -H° K ). 



(20) 



To coarse-grain Hk, we use the following projection 
operators, 

C(A7r'-) = h i0 n^...n T jt a jt^...(l - n ir ), (21) 

where the Greek letters refer to (a, a), where a = (1,2) 
refers to the orbital degrees of freedom, and a ^ ^ 
7'''- P JaWi 7> ' ' ' ) P ro j ec ts out wavefunctions which 
have n electrons occupying states ((3, 7, • • • ) and no elec- 
trons in the remaining states on site i except a. With the 
help of these projection operators, we can rewrite Hk as 



2n a 



H K 



E**E E E 

(ij)aba n,m=0 (/3/y,...) (#',7',— ) 

^" 0lCT )(/3,7, • • • )4a^?M)(/^7', • • • )■ 



(22) 

The hnal step is to classify the terms in Eq. |2"2"1 according 
to the change of interaction energy Ej ^ 0. This can be 
determined by the commutation relation, 



[HK,E r , Hi] — —EiHk,E!- 



(23) 



Substituting Eqs. [23] into [2(2 we can express Si in the 
general form 



Si = 



H 



Ei 



K,E, 
E, 



(24) 



Although the above formalism is quite general, it is 
cubersome to write out all the terms explicitly for n - 
orbital systems. In particular, if we are interested in the 
low energy effective theory, only parts of Si are neces- 
sary depending on the filling. In the next section, we will 
present the results for first-order corrections to the inten- 
sity of the lower energy bands, m h (n), for filling n = 1—x 
as an example. 



B. Intensity of lower energy bands for n — 1 — x in a 
two-orbital Mott system 

Terms that must be retained in Si for n = 1 — x are 
those associated with hopping processes that create or 
destroy doubly occupied sites, namely 



h' k = - >: + ^ 



E 

(ij)ab<T /3 



(25) 



The evaluation of the commutator [H' K ,Hi] is compli- 
cated by Hund's coupling and the paired hopping. The 



recipe for treating the case with a non-zero J is as follows. 
We consider the term 



*i-f/(if) [(1 i)] Ciit^iif-P/(it) 



(26) 



which describes the hopping process of \il t) ~~ ^ 

\il I 1 t,j0). Because \il l 1 f) is not an eigenstate of 
Hi (see Table [J) , a commutator of the form in Eq. [23] 
can not be obtained. Instead, we need to introduce 



A-ij = P i(if) [(1 i)]%t 6 jit-Pj(it) 

±4^H^W) [C 1 D] 4 t ^'it^(it) • ( 27 ) 
Straightforward calculations give 

[At j ,Hi}\ip) = {At j H I -HiAt j )\i>) 

= (4# z -#z4)|iU,ilt)<8|...> 

= -ff / [|iUit,io)±|i2i2t ! io)]<8|...) 

= -([/±J)[|iUlt,jO)±|i2i2t,jO)] ®|...) 



-(U±J)AfM). 



(28) 



In the above equations, we have used the fact that Aij 
only changes states related to sites i and j in l^) so that 
we just need to focus on the change of interaction energy 
on sites i and j. Combining Eqs. [20l and [28l we find 



-U 



A+ 



A; 



<U + J U-J 



),H I ]=t 1 Pl {lt) [(li)]c} lt c m P^ y 

(29) 

Following the same procedure, we organize all the terms 
in the first-order correction associated with H' K as 



S[ 



Sl - E i ~ gV^gj 

1=1 



Ei 



(30) 



where E 1 = U'-J, E 2 = U'+J, E 3 = U-J, E 4 = U+J, 
and si_-e, = s \ Er Detailed expressions for si^e, are 
given in Appendix A. 

With these results in hand, we are now ready to evalu- 
ate the intensity of the lower energy bands. The unitary 
transformation 



4aa + [ S ^4aa] 



(31) 



relates the bare electrons to the transformed electron op- 
erators. To derive the spectral properties, we need to 
separate Eq. [3T]into operators which are coarse-grained 
relative to the energy scale. In particular, the zero energy 
component yields 



J = gt 



4 1 

i=i 1 



where 

°iatr,0 
c iacr,U'-J 



a f p° 

n ' "" " ~^ ^iba^ibS^iaS , (32) 
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and b ^ a. Using the approximate expression for the 
intensity of the lower energy bands derived from the first 
moment of the spectral functions 1 - , we have 

m h (l-x) = -^^2({ciaa,o,cl aafi }) 

iaa 
iaa 

4 1 

~ ^E (-^{lc™°,E n ,si, En ],c\ aafl } + h.c.)) - 



While the first term in Eq. [33] yields 1 + 3x which is 
simply the spectral weight transfered in the atomic limit 
given by Eq. [9] the second term in Eq. [33] is precisely 
the DSWT in question. This expression can be simplified 
using two relationships. First, Ci a< j,E\G) = {G\c\ aa E = 
since we are considering the case with filling n = 1 — x < 



1. Second, c iaa fiCia<j,E 
these properties to Eq. [I 



: C-ia<7.0c\ a 

I leads to 



0. Applying 



... . {{[^aa,E n ,S^ En ],c\ a(j0 }) = {Ci aa , E Sl, E c\ a<rfi ), (34) 

(3ai)d a straightforward calculation yields 



(Cia<7,U'-jSl,U'-jctaa,o) = ^M^ibtr^+SM) + ^{c^h+S^a) + E ^ (4><r ^+S,ba) + "TjT (4ba S i+6,air ) > 

5 5 

-(Cia<T,U<+jSl,U>+jCia<T,o) = E ^ (4b9^+S,bS ) + "y (c\b9 £ i+S,aa) , 

S 

i~ ~t \ * a /~t - \ i /-t - \ 

-\Ciaa,U-jSl,U-jCia<T,o) - / , TT \ c iaS C i+SMS} + ~ \ c ia<j C i+S,ba) ' 

-(CtaaM+jSl.U+jclaafi) = ^laS^+SMa) + ^ {c\ aS h+S,ba) ( 35 ) 



where 5 = x,y. Substituting Eq. 
finally arrive at 



m h (l - x) 



1 + 3x + a, 
_/ 1 



into Eq. [33J we single-band Hubbard model to 2 ( 



U'-J 



u' 

U r2 -J 2 



u - 

' U 2 -.] 2 



U' 



U 



a 



J U' 2 -J 2 U 2 -J 2 



(36) 



where a is a measure of the DSWT with 



f 



M 5 

<l,J>,<7 

17$ 



(37) 



the average kinetic energy of the empty site created in 
the lower band per site evaluated in the ground state. 
Note the presence of transformed operators in the ex- 
pression for Hk ■ This result represents a generalization 
of the single-band Hubbard result^^ to the two-band 
case in the presence of Hund's coupling. It can be seen 
that DSWT is enahnced in the two-orbital system be- 
cause the prefactor of Kk in a shifts from 2/U in the 



The physics behind this trend is as follows. Since DSWT 
is induced at the expense of creating doubly occupied 
sites, more orbitals provide more ways to create such 
sites, which results in the enhancement of DSWT in a 
multi-orbital case. 
C. Effective low energy Hamiltonian for n = 1 — x 

While it is certainly complicated, we can also derive 
the effective Hamiltonian in the transformed basis. Eq. 
[17] leads to the effective Hamiltonian 



H = H K + H!+ [Si,H K ] 



(38) 



Given the explicit expression for Si for n = 1 — x, one 
readily finds that 



eff 



H 



where 



E 



b^a 2,(5,5' ,<r 



-IT 1 



Ha 



J U' 2 - J 2 U 2 - J 2 

(39) 



-h.c.) 
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Hi = -t 2 a Pi+S' ,(aa)Pi+8,(aa)^l+S' ,aa^+S,aaPi^ a(T )[(bcr)] - tatbPi + S' : (bcr)^i+S' ,bcr^ ibaP l(a<y) l(b&)]c\ a(7 Ci + 8Mcr Pi + S,(aa) 

H J 2 = -^H^a^+<5^(aff)^^+^,(a(T)cl+^^a(TC^+5,a(T-f^ 1 ( a(T )[(^0•)] + i a ifcP" +(5 , | ( fcs ) c\ +s , bs CibS P^aa) [(b^)}cl aa C t+ S,acrPi + S,(aa)} 

~ J^aPi+S' ,(acr)Ci+S' ,aa^ i , a ^^i,b<7^ i - bS ^i^0^ i + S ^ a < T Pi+S,(aa) 

- Jt a tbPi +g , t (ba)Ci+8' ^,6CT^CTC i+5)6CT Q+5,bCT^/ + 5 ) ( Q CT)[(^)]Cj + 5 ia CT^aCT^°( aCT ) 

^3 = _ ^{ i a^°+5\(aa)^+<5 : (a ( T)Ci + 5^ a o'^+^aff- P il(acr)[( a ^)] + *a-^+5',(aff)^!+5',aCT^aCT f£( ocr ) [(^)]cJa<r^H-i,aff-^4.5,(a<r) } 

— 7> 2 P° r- 7$ o. -7& r- r P° 

~ Jta\Pi+8' ,(a5)c\+8' ,aaCiasc\ + $ £ S C~i+&,aoPi+S ,(ba)[( a ^')}^i+S,b(T^ ia ' 7 ^( aa ) 

^4 = ^\2Pl+8,(ba)\( aa )]Pi+8' ,{aa)^i+S\aa^i+^,aaPi + s^aa) ~ *a*12-Pj°|_,5' , ( a(T ) Cj +(J , , acr Cib<j Pl^ a<J ) [(b&)]cl aiJ Ci + 8,aa P° + s ,( a(7 ) 

#5 = _ ^'*12-Pj+5,(6ct) [( a ^)]^i+<5',(acr)c! + 5' i a ( TC l +5,a ( TP° + 5,( a(T ) 

#6 = ^^ ^ L^^+<5,(6 CT )[(^)]^+<5^(a CT )ci + 5' IaCT C ^+ ^ Ia(T P i t _ (5J(Qcr ) 

(40) 



In deriving this expression, we used the properties 
c W-Fi( a(7 )[(fc r ')]|£r) = and Hi\G) = since the ground 
state can only have singly-occupied sites for n = 1 — x. 
In Eq. [40l three-site hopping terms can be identified 
with 6^6'. Moreover, terms with 6 — 6' represent 
generalized superexchange interactions between nearest- 
neighbor sites occuring not only between spin but also 
orbital degrees of freedom. Several terms proportional to 
J are induced as a result of the internal change of states 
on the same sites allowed by the Hund's coupling and 
paired hopping. As was the case with the single-band 
Hubbard model, this effective model contains three-site 
terms which can not be obtained by simply considering 
generalized superexchange processes between spin and 
orbital degrees of freedom. 

D. Particle- hole symmetric Quantum critical point 

Recently, it was proposed 1 ^ that the thermopower of- 
fers an unbiased scale for determining the doping level in 
the cuprates. The basis for this observation is that the 
thermopower for a wide variety of hole-doped cuprates 
all collapse onto a single curve and crosses zeroii at a 
universal doping level of x = 0.23. Since there is only 
one fitting parameter in the thermpower and it is the 
doping, the thermopower can be used to map oufeii the 
high-T c phase diagram. Physically, the thermo-electric 
power S is defined as the voltage drop across a material 
generated by a temperature gradient. Since itinerant car- 
riers flow from higher to lower temperature regardless of 
their charge, the vanishing of S occurs when particle-hole 
symmetry is present. 

In Fermi liquid theory, particle-hole symmetry only 



holds at half-filling and any finite doping of holes should 
break this symmetry. As a result, the explanation for a 
vanishing of S at optimal doping requires the incorpora- 
tion of strong correlation effects. One of u s 18 i 19 has ad- 
dresed this issue by considering the spectral weight trans- 
fer in the single-band Hubbard model in the strong cou- 
pling limit. The particle-hole symmetric condition can be 
defined where the number of states below and above the 
chemical potential becomes equal. In the atomic limit, 
this condition yields 

l-x c = 2x c , (41) 

which predicts the vanishing of S at x c = 1/3. Away 
from the atomic limit, DSWT renders x c to a lower value 
around x c = 0.23 as observed in the experiments. That 
particle-hole symmetry is dynamically generated near op- 
timal in strong coupling models of the cuprates has also 
been confirmed by Vidhyadhiraja and co-workers 20 . 

This senario can be generalized to the current n - 
orbital Mott system, allowing us to make the following 
predictions. We present the particle-hole symmetric con- 
dition in the atomic limit to show how the critical doping 
with vanishing S depends on the multi-orbital structure, 
and it is expected that in the real materials this critical 
doping should be reduced to lower values due to DSWT. 

In the case of n = 1 — x, the particle-hole symmetric 
condition is unambiguously given by 

1 - x c = 2n a x c =£• x c — - — - — - (42) 
2n D + 1 

in the atomic limit which is strictly less than 1 /5 as would 
be the critical doping in the 2-orbital case. Interestingly, 
Mukerjee^i has analyzed the thermopower in the atomic 
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limit of the 3-orbital Hubbard model and found it to 
vanish at the doping level predicted by Eq. (|4"2"j) . This 
agreement further corroborates our treatment here that 
the particle-hole symmetric condition in the lower band 
shifts to lower values of doping as multi-orbital contri- 
butions are included. For n > 1, again the situation is 
complicaed by Hund's coupling and the paired hopping, 
and the particle-hole symmetric condition depends cru- 
cially on the magnitude of J. If J is large enough to 
produce new gaps in addition to those created by U, the 
particle-hole symmetric condition should be determined 
only by the bands in the vicinity of the chemical poten- 
tial. For example, for the case of n = 1 + x shown in Fig. 
(He), the particle- hole symmetric condition should be 

2x c = 1.5(1 - x c ) => x c = j (43) 

if J is large enough. If J is not large enough to generate 
new gaps, all the states within fi — E~ and fi + E + defined 
in Eqs. [TT] and [T^] should contribute to the thermopower. 
Consequently, in the latter case the particle-hole sym- 
metric conditions can be written in general form 

n*(l -x) = (2n - m + l)x => x h c {m) = tt^ 1 — (44) 

2n + 1 

for n = rii — x, and 

2tI — Tl ' 

{n i + l)x = {2n -n i ){l-x)^x%{n i ) = ° / (45) 

2n a + 1 

for n — Ui + x. 

These relationships have interesting consequences for 
superconductivity. Since the paricle-hole symmetric crit- 
ical doping is intimately related to the optimal doping of 
superconductivity^—, we can draw the following con- 
clusion. Optimal doping for superconductivity in multi- 
orbital Mott systems should be reduced relative to their 
single-orbital (that is, one unpaired d electron) counter- 
parts because the particle-hole symmetric conditions de- 
rived above are all bounded from above by the single- 
orbital value of x = 1/3. According to Eqs. l44l and 
@5J cither a hole-doped multi-orbital Mott material with 
smaller integer fillings, or an electron-doped Mott multi- 
orbital material with larger integer fillings could serve as 
a good parent compound for new high-T c superconduc- 
tors. 



weight of the unnoccupied part of the lower band grows 
as 4x as opposed to 2x in the single-orbital case. This re- 
sult can be easily be confirmed by doping a multi-orbital 
Mott system such as NiO and measuring the oxygen K- 
edge x-ray intensity as has been done previously for the 
cuprates 3 . As the hopping terms are turned on, the re- 
distribution of the spectral weight, namely dynamical 
spectral weight transfer (DSWT), occurs at all the en- 
ergy scales in a manner similar to that in the single-band 
Hubbard model. 

We have developed a general formalism to perform the 
perturbation theory in the strong coupling limit and ap- 
plied this method to the two-orbital system. The inten- 
sity of the lower energy bands has been explicitly calcu- 
lated for filling n = 1 — x, and we have found that the 
leading-order correction due to DSWT is proportional to 
the average kinetic energy of the empty sites. Relative 
to the single-band Hubbard model, we have shown that 
both static and dynamical spectral weight transfers are 
greatly enhanced as a consequence of multi-orbital struc- 
ture. The corresponding low energy effective Hamilto- 
nian for n = 1 — x in the two-orbital Mott system has 
been derived. Moreover, we have shown that a critical 
doping, defined as the doping level exhibiting zero ther- 
mopower, is reduced to a smaller value as the number 
of orbitals increases. Motivated by the lesson from the 
cuprates that this critical doping is intimately related to 
the optimal doping of superconductivity, we have pro- 
posed that either a hole-doped multi-orbital Mott ma- 
terial with smaller integer fillings, or an electron-doped 
Mott multi-orbital material with larger integer fillings 
could serve as a good parent compound for new high-T c 
superconductors. 
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VI. APPENDIX A 



IV. CONCLUSIONS 

In this paper, spectral properties and related predic- 
tions have been studied in multi-orbital Mott systems 
for the first time. Taking a two-orbital system with full 
on-site interaction as a starting point, we have analyzed 
exact spectral functions in the atomic limit near inte- 
ger fillings of n.i = 1,2,3 and generalized the results to 
n Q orbitals at arbitrary filling. Most striking is the re- 
sult that for n = 1 — x, the static part of the spectral 



In this appendix, we write out S[ explicitly following 
the same procedure described in the Section IIIIB1 To 
simplify the notation, we use a short-hand index of 1 — > 
(l,t), 2 -> (1,|), 3 -> (2,t), 4 -> (2,|). This leads to 
the following result: 

4 = pA(2)4%i^±4«Jf 3 (2)42;if?i. 
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Br 



' ■ ),Hj] =-^(1)4^2, 



Lo 



2 y U + J U-J' 



C 



± _ 



pUm^p% ± (4)4 c j3 pf 3 



C-. 



2M + J U-J 



Pu{3)Ci4pj4,Pj4 ^ 4^3-Pi2(3)c i2 Cj4-Pj4) 

d: 



L 2 V P + J ^ U-J' 



Mi] =-PIMcU 3 aP% 



E% = ^ 1 3 (4)4c jl P? 1 ± 4^(4)4 c iX Pl 



e: 



2M + J U-J 



),H I \ = -pum^iPh 



PlMc\^2P% ± 4^3^(3)4^2^2, 

F " \Hi]=-p!a(M^p% 



1 P + 



G, 



2M + J U-J' 

PA(2)45 J3 ^° ± clc i2 P? 3 (2)cl 3 c j3 P? 3 , 
G 



i 1 ( G ti 

[ 2 { U + J ^ U - J' 



,H I \=-Pk{2)c\ l c j3 P%, 



Pai^^CjiPji ^ ^i 3 ^Pii(^)^i4^jiPj4y 

[1 ( B$ , Hi 



2M + J U-J' 



Mi] = -P i 1 2 (l)4^4P°4- 



We can also apply the same technique to hopping terms 
leading to final states with \il t 2 i,jO) or \il J, 2 f 
, jO), but this time there exists two different ways. For 
example, we can define, 



with t±2 in this category without loss of generality. We 
find immediately that 



K 



pA (4)4^i^i ± 4^a (4)4^-1^1, 

Lo V TV 



2M'-J U' + J 



)Mi] = ^(4)6^^°, 



P^ml^P^ ± 4^(3)4^. 



■1, L 



-2 y U'-J U' + J' 



• H I ]=-Pl 2 (3)~cj 2 d j2 P^ 



M± = PA(2)4c, 3 P^± 4^(2)4^, 



■1, ilp 



AP 



-2 y U'-J U' + J' 



Mr 



-Ph^)c\ 3 c 33 P% 



N 



± _ 



^AW^^P^ ± C l2 C l iP l 1 3 (l)Cj 3 C J -4Pj ) 4, 



2 K U'-J U' + J'' 



Hi] = -Pk{l)c\ A c 3i P% 



0% = P i3 (2)45 jl P°T 4^2^(2)4^^, 



-I, o 



0; 



-2 K U' - J U' + J' 



Mi]=-Pl 3 {2)c\ 3 c jl P^ 



Q% = P/4(l)4^2P°2T4^iA 1 2(l)4s j 2Pj 2 , 



-1 



[o \TT> 



2 U' — J U' + J 



)Mi] = -pLm^ 2 p° 2 , 



R 



± _ 



Pam^aPis t4^«(4)45j3^3. 

r 1 / P^ Pi 



-2 y U'-J U' + J' 



Mr 



^(4)4^3. 



P± = PA (4)4^-iP° ± 4^4^2(4)4^1^1, (46) VK± ee PA (3)4^4^ T 4^(3)4^4. 



K% = PA(4)4^ifi T 4^^3(4)4^1^1. (47) 



L 2 V P'-J U' + J' 1 



Hi] = -P/ 2 (3)ct 2 c^P 



J'4- 



Both choices can lead to, 



■1, ^ 



K: 



L 2 y U'-J U' + J' 



Mi] = -P i \(4)cJ 1 c jl P° , (48) 



This is because for the transformation from |1 t 2 J,) to 
1 1 I 2 t) , we can either flip the spin (Eq. or exchange 
the orbital (Eq. |47|) . Since both choices are legitimate 
and related up to another unitary transformation, we can 
choose the spin-flip scenario for the hoppings with t\ and 
t 2 and the orbital exchange scenario for the hoppings 



For hopping terms leading to final states with \Ha2a,jO), 
no complication occurs because these final states are still 
eigenstates of Hi. It is then straightforward to obtain 

Tv-j= E - M^i (3)4^i^i +^2(4)4^2^°] 

<i,j> 

- *2 [Pl 3 {l)c\ 3 C l3 P% + PkV)4&*P%] 

- - tia [PUM&iPjx + Pkm~clc 32 P% 

+ P i 1 1 (3)4c J -3Pj 3 + PA(4)4^4P, C 4 ] • 

(49) 
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Collecting all the above terms, for n/ < 1, we find that where E\ = U' — J, E 2 = U'+J, E3 = U — J, E4 = U+J, 
the relevant terms in S-t is and s-i _ w — s\ ^ . 



S'x - E 



4 

Si E, - S\ —Ei 



1=1 



E, 



81,V-J = f W -j+ E -*l(^ + i$)-*2(M± + JV±) 
<i,j> 

<i,j> 

= E -*Mi + B«)-* 2 (C^ + !>,-•) 

<j,j> 

t 12 (S i T + F i T + G- +J ff-) 

si.u+j = E + B tii - HCtj + ^5) 

<i,j> 

(50) 
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